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The ratios of charged antiparticles to particles have been 
obtained for pions, kaons, and protons near midrapidity 
in central Au+Au collisions at ^s NN — 200 GeV. Ra- 
tios of <7r-)/(7r+) = 1.025 ± 0.006(stat.) ± 0.018(syst.) , 
(K-}/(K+) = 0.95 ± 0.03(stat.) ± 0.03(syst.), and <p)/<p) = 
0.73 ± 0.02(stat.) ± 0.03(syst.) have been observed. The 
(K~) I (K + ) and (p)/(p) ratios are consistent with a bary- 
ochemical potential /j,b of 27 MeV, roughly a factor of 2 
smaller than in ^/s NN = 130 GeV collisions. The data 
are compared to results from lower energies and model 
calculations. Our accurate measurements of the particle 
ratios impose stringent constraints on current and future 
models dealing with baryon production and transport. 

PACS numbers: 25.75.-q 

In this paper, the ratios of antiparticles to particles 
for primary charged pions, kaons, and protons in central 
Au+Au collisions at y/s NN — 200 GeV were determined 
using data from the PHOBOS detector. The data were 
taken during the 2001 run of the Relativistic Heavy- 
Ion Collider (RHIC) at Brookhaven National Laboratory. 
The experiments at the RHIC aim at understanding 
the behavior of strongly interacting matter at high 
temperature and density, testing predictions of quantum 
chromo dynamics (QCD). 

Earlier studies showed that the (p) / (p) ratio increases 
from 0.1 observed in fixed-target experiments at ^/s NN = 
17.3 GeV [1,2] at the SPStoavalue of 0.60±0.04(stat.)± 
0.06(syst.) observed in PHOBOS at ^/T^ = 130 GeV 
[3] which is in agreement with the values observed by 
the other RHIC experiments [4-6]. Nucleus + nucleus 
(A+A) collisions at the RHIC have reached a regime 
where the fraction of baryons at midrapidity is not 
dominated by baryon transport alone. Therefore, these 
collisions are approaching the net baryon-free regime for 
which experimental results can be directly compared to 



first-principle QCD calculations. The particle ratios at 
130 GeV were shown to be consistent with statistical 
models of particle production [7], assuming chemical 
equilibrium. These models gave a very good description 
of the particle ratios with an assumed hadronization 
temperature of 175 MeV and a baryochemical potential 
He =45 + 5 MeV. However, it proved difficult to explain 
the results in terms of microscopic descriptions of the 
evolution of the particle source using hadronic transport 
models [3,8,9]. 

The data in this paper allow further tests of models 
of baryon production and baryon number transport in 
nuclear collisions, using a study of the energy dependence 
of particle ratios between y/s NN = 130 GeV and ^s NN = 
200 GeV. The results were obtained using the PHOBOS 
two-arm magnetic spectrometer. Details of the setup 
have been previously described [10-12]. Each arm has 
a total of 16 layers of silicon sensors (corresponding 
to 50000 pixels), providing tracking both outside and 
inside the 2-T field of the PHOBOS magnet. For this 
analysis, only the central region of each arm was used, 
covering ±15° around the 45° axis of the experiment 
in the horizontal plane. As three silicon layers do not 
extend into this geometrical acceptance region, particles 
within this acceptance region typically traversed 13 sili- 
con layers. A two layer silicon detector (VTX) covering 
\rj\ < 1.5 and 25% of the azimuthal angle provided 
additional information on the position of the primary 
collision vertex. The two arms and frequent magnetic 
field polarity reversals allowed for many independent 
checks in the determination of particle ratios. The use of 
highly segmented silicon detectors, beginning at 10 cm 
from the interaction point, result in systematic errors in 
the correction to the ratios from secondaries and hyperon 
decays of 2% or less. 



The primary event trigger was provided by two sets of 
16 scintillator paddle counters, which covered a pseudo- 
rapidity range 3 < \rj\ < 4.5. Additional information 
for event selection was obtained from two zero-degree 
calorimeters which measured spectator neutrons [13]. 
Details of the event selection and centrality determina- 
tion can be found elsewhere [14,15]. Monte Carlo (MC) 
simulations of the apparatus were based on the HIJING 
event generator [8] and the GEANT 3.21 simulation 
package, folding in the signal response for scintillator 
counters and silicon sensors. 

For this analysis, events with the 12% highest cross 
section were selected. A Glauber calculation was used 
to relate the fraction of the cross section to the number 
of participating nucleons, N part [14,15]. The average 
number of participants for these events was estimated 
to be (Np ar t) — 317 ± 10(syst.). The data used in this 
analysis was taken during the first two months of RHIC 
2001 running and represents only 10% of the total data 
set. 

The geometrical layout of the PHOBOS detector led to 
an asymmetry in the acceptance and detection efficiency 
for positively and negatively charged particles, h + and 
h~, at a given magnet polarity, B+ or (see Fig. 1). 
Particle and antiparticle yields obtained with opposite 
field polarities were used to calculate the (tt~) / (tt + ), 
(K-)/(K+), and (p)/(p) ratios as in Eqs. (1). 
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It is assumed in the analysis (and checks were performed 
to verify this) that the acceptance and efficiency for a 
specific particle type at P> + is identical to the acceptance 
and reconstruction efficiency for the antiparticlcs of that 
type at B _ , and vice versa. This measurement requires 
the further assumption that the pr dependence of the 
particle and antiparticle yields within our acceptance is 
similar. 

The reproducibility of the absolute field strength was 
found to be better than 1%, based on Hall probe mea- 
surements for each polarity and the comparison of mass 
distributions for identified particles for the two polarities. 
The trigger event selection was checked using the number 
of reconstructed straight line tracks in the low magnetic 
field region. These numbers agreed to within 0.35% for 
the B + and B~ data sets. An asymmetry in either the 
total field strength or the trigger selection would shift one 
of the measurements of each ratio in Eqs. (1) up while 
the other would shift down. Since each shift is equal (and 
small), this divergence between the two measurements for 
each particle type in Eqs. (1) cancels in the average. 



Only events with a reconstructed primary vertex posi- 
tion between -6 cm < z vtx < 6 cm along the beam axis 
were included in this analysis. This vertex selection en- 
sures that both particles and antiparticles can be tracked 
and identified in the spectrometer for both bending 
directions, and thus reduces vertex-position-dependent 
systematic effects. The vertex position was determined 
by two methods. One method combined information 
from both spectrometer arms and the vertex detector 
("global vertex") and the other only used information 
from the arm for which the ratios were being measured 
("local vertex"). The global vertex is more precise, but 
can be affected by misalignments between the detectors 
while the local vertex is less precise, but insensitive to 
misalignment between the different detectors. 
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FIG. 1. Contours of the acceptance of the spectrometer as 
a function of transverse momentum and rapidity for pions, 
kaons, and protons, where the acceptance has fallen to 10% 
of its maximal value. The left plot is for particles bending 
towards the beampipe (h + ,B + or h~,B~) and the right plot 
is for particles bending away from the beampipe (h + ,B~ or 
h~,B + ). The acceptance is averaged over the selected vertex 
range and the accepted azimuthal angle (f>. This results in pr 
acceptance ranges of 0.1-0.6 GeV/c for 7r +, ~, 0.1-0.5 GeV/c 
for K + '~ and 0.15-0.9 GeV/c for p,p. 

The steps of the track reconstruction algorithm have 
been previously described [3]. For the present analysis, 
an additional track fitting step was added, which takes 
into account the correlations between hit positions along 
the track. The final track selection was based on the x 2 
probability of this fit, improving the rejection of tracks 
with incorrectly assigned hits and thereby improving 
momentum and particle identification resolution. 

Particle identification was based on the truncated 
mean of the specific ionization dE/dx observed in the 
silicon detectors. The identification cuts for pions, kaons, 
and protons are shown in Fig. 2. The corresponding 
acceptance regions for identified particles in transverse 
momentum px and rapidity are shown in Fig. 1. 

The accuracy of the ratio determination relies on the 
precise cancellation of acceptance, efficiency, and back- 
ground effects between the two magnet polarity settings. 
To allow a detailed cross-check of this cancellation, the 
data were divided into statistically independent subsets 



(see Table I). The first division was necessitated by a 
significant transverse shift in the beam orbit (0.5 mm) at 
one point during data taking. Data before and after this 
shift were analyzed separately. The data from the two 
spectrometer arms were also analyzed separately. This 
resulted in four statistically independent subsets. These 
four subsets resulted in eight independent measurements 
of each ratio (each bending direction in the magnet gives 
a separate measurement). Each data set was analyzed 
using both the local and the global vertex position, 
resulting in a total of 16 measurements of each ratio. 
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FIG. 2. Distribution of average truncated energy loss as a 
function of reconstructed particle momentum. Three clear 
bands can be seen, corresponding to pions, kaons, and 
protons. The solid lines indicate the cut regions for counting 
identified particles. 

Table I gives the event and identified particle statistics 
for all combinations of magnet polarity, particle charge, 
particle species, spectrometer arm and beam orbit con- 
dition. The table contains the results using both the 
global vertex and the local vertex as the starting point 
for the track finding. The systematic uncertainty due 
to the vertex finding method was extracted from the 
distribution of half of the difference between the value 
found using the global vertex and the value using the 
local vertex. The rms of this distribution, used as an 
estimate of the systematic uncertainty, was found to 
be ±0.003 for ±0.018 for {K~)/{K + ), and 

±0.013 for (p)/(p). 

To perform consistency checks, Eqs. (1) were used to 
calculate (7r~)/(7r + ), (K~)/(K + ), and (p) / (p) for each of 
the eight statistically independent data sets. To check if 
a systematic deviation exists between ratios calculated 
in separate spectrometer arms, a weighted average of 
the four ratios calculated in the positive spectrometer 
arm was compared against a weighted average of the 
four ratios calculated in the negative spectrometer arm. 
Half of this difference was taken as an estimate of 
the systematic error associated with ratios calculated 
in separate spectrometer arms. As the pion ratios 
have the smallest statistical error, they possess the best 
ability to resolve systematic deviations. Therefore, if 
this difference for kaons or protons does not exceed a 
one sigma statistical deviation, half the pion difference 



was used. This same method was performed to estimate 
systematic errors associated with the other two divisions 
of the data. The results of these checks showed that no 
deviation was evident beyond the statistical errors for 
kaons or protons so the pion ratio's estimate of ±0.015 
for the beam orbit division, ±0.001 for the arm division, 
and less than ±0.001 for the bending direction division 
was applied to the kaons and protons. 

Table I also shows the average transverse momentum, 
(Pt), of identified particles, not corrected for acceptance 
and efficiency. For all particle species, the (px) of 
accepted particles at one magnet polarity agree to within 
±3% with the (pr) of the accepted antiparticles at the 
opposite polarity. Also, the average rapidity and the 
second moments of px and y distributions for accepted 
particles and antiparticles agree within errors. As a 
consequence, no acceptance corrections were made to the 
particle ratios. 

TABLE I. Number of accepted events, accepted identi- 
fied particles and average transverse momentum ({pt}) in 
(MeV/c) of these particles at each magnetic field polarity. 
Data are split up by beam orbit condition and spectrometer 
arm. The number in parentheses is the number of particles 
found using the local vertex, while the number not in 
parentheses uses the global vertex. 



Beam orbit condition 1 




B+ 


23725 events 


B" 


■20965 events 


Narm 


-i* particles 


(Pt) 


N particles 


(Pt) 


7T + 


8532 


(8286) 


418 ± 1 


12501 


(12935) 


289 ± 1 




14298 


(14692) 


288 ± 1 


7643 


(7472) 


416 ± 1 


R+ 


200 


(195) 


379 ±4 


241 


(255) 


259 ±4 


R~ 


243 


(250) 


253 ±4 


208 


(198) 


378 ±3 


V 


458 


(472) 


580 ±6 


396 


(405) 


481 ± 6 


V 


319 


(324) 


480 ± 7 


293 


(291) 


596 ±8 


Parm 


^ particles 


(jpr) 


^ particles 


(Pt) 


7T + 


7139 


(7273) 


414 ± 1 


13172 


(13423) 


279 ± 1 




14961 


(15344) 


278 ± 1 


6432 


(6569) 


413 ± 1 


R+ 


215 


(223) 


373 ±3 


280 


(290) 


256 ±3 


R~ 


308 


(310) 


256 ±3 


179 


(182) 


373 ±3 


P 


409 


(427) 


584 ±6 


424 


(433) 


485 ±6 


P 


327 


(346) 


474 ± 7 


267 


(264) 


573 ± 7 


Beam orbit condition 2 




B+ 


28785 events 


B~ 


: 13824 events 


Narm 


N particles 


(Pt) 


-i* particles 


(Pt) 


7T + 


10471 


(10248) 


417 ± 1 


7407 


(7759) 


290 ± 1 




16607 


(17267) 


291 ± 1 


5083 


(4990) 


417 ± 1 


K+ 


304 


(302) 


376 ±3 


137 


(154) 


257 ±5 


R- 


278 


(305) 


257 ±3 


135 


(124) 


373 ±4 


P 


612 


(606) 


577 ±5 


248 


(262) 


492 ±8 


P 


398 


(407) 


484 ± 6 


195 


(193) 


573 ±9 


Parm 


N particles 


(Pt) 


^ particles 


(pt) 


7T + 


8322 


(8574) 


414 ± 1 


8746 


(8946) 


280 ± 1 




18680 


(19181) 


276 ± 1 


4229 


(4350) 


413 ± 1 


K+ 


238 


(249) 


369 ±3 


202 


(193) 


260 ±4 


R~ 


396 


(398) 


257 ±3 


106 


(104) 


378 ±4 


P 


542 


(565) 


584 ±5 


289 


(298) 


479 ± 7 


P 


418 


(434) 


491 ± 6 


179 


(197) 


584 ±9 



After the systematic uncertainties of the raw particle 
ratios were determined, the weighted average and sta- 
tistical error for each of the three ratios was obtained. 
The errors associated with vertex finding method and 
consistency checks were added in quadrature to give an 
estimate of the total systematic error. The result was 
±0.015, ±0.02, and ±0.02 for pions, kaons, and protons, 
respectively. 

Contamination of the proton and kaon samples by 
misidentified particles was checked by testing the stabil- 
ity of the particle ratios against variations of the particle- 
ID cuts shown in Fig. 2. Within the statistical error the 
ratios were stable against changes of these cuts. The 
effect of pion contamination on (p) / (p) was estimated by 
modeling the pion's truncated energy loss distribution 
as a function of momentum. A conservative estimate 
of the number of pions falling within the proton and 
antiproton cut region produced a correction of less than 
0.8% to (p)/(p), which is reflected in its systematic error. 
This same correction for (K~) / (K + ) was estimated to 
be less than 0.1%. For pions, electron contamination 
to (ir~)/(n + ) was estimated to be less than 0.2% based 
on studies carried out using GEANT simulations with 
HIJING events. 

To represent the ratios of primary particles produced in 
the collision, the detected particle ratios have to be cor- 
rected for particles produced in secondary interactions, 
loss of particles due to absorption in detector material, 
and feeddown particles from weak decays. Secondary 
and feeddown corrections were described in detail in the 
y/s NN — 130 GeV paper [3]. As in that analysis, these 
corrections were minimized by requiring tracks to point 
back to the vertex within 3.5 mm. The total effect of 
all corrections to the pion and kaon ratios was estimated 
to be less than 1%; this correction is reflected in the 
systematic error. 

Corrections to (p)/(p) are more significant. The 
secondary correction from the production of protons in 
the beampipe and detector material was found to be 
0.7%. This correction is very small in the spectrometer 
for the following reasons. Secondary protons tend to be 
produced at low momentum. High multiple scattering 
at low momentum makes detection efficiency of such 
secondaries small. Particles are also required to pass 
through at least 11 silicon layers. Low momentum 
particles tend to bend out of the acceptance (as can 
be seen from Table I, where the (pr) of protons and 
antiprotons is above 450 MeV) . The number of secondary 
particles produced scales with the number of primary 
particles produced, mostly pions. The correction is 
calculated using HIJING events and, therefore, depends 
on how well HIJING reproduces the measured p/ir ratio 
and their momentum distributions. Measurements of 
p/ir at 130 GeV [5] are consistent with the values from 
HIJING at the same energy. A systematic error of the 
size of the correction (1%) is conservatively assigned. 

The absorption correction arises from an asymmetry in 
the loss of antiprotons versus protons interacting in the 



beampipe and planes of the spectrometer. A momentum 
dependent correction is used. Only GEANT's ability 
to model interactions in material plays a role in this 
correction. The hadronic interaction packages Gheisha 
and Fluka were both tested. Gheisha produced a 
correction of 5.1% while Fluka produced a correction of 
2.2% The correction was therefore estimated from the 
average as 3.7% with a systematic error contribution of 
1.6% taking into account the errors on each individual 
package's result and the difference between those results. 
This correction was found to be consistent with analytic 
calculations. 

The feeddown correction is calculated from Eq. (2): 

(p\ /( p+p^A = 1+ ^-t(f/f) 

\p)/ \P + P A ^-J l + br-e acc -± 

This accounts for the difference in the number of hy- 
perons (primarily A and A) decaying to protons versus 
antiprotons. e acc is the fraction of protons from lambdas 
or antiprotons from antilambdas detected within our 
acceptance and br is the branching ratio for A — > pir~ 
or A — > pTr + . br ■ Eacc equals 0.3 averaged over px- The 
Pt dependence of e acc is weak due to two competing 
effects at high and low momentum. Higher momentum 
hyperons decay later due to time dilation. The decay 
particles from these hyperons have a reduced chance of 
being reconstructed because they are more likely to decay 
past the first planes of silicon in our spectrometer. Lower 
momentum hyperons have decay particles which deviate 
more from the original path of their parent and, therefore, 
have a lower chance of pointing back to the primary 
vertex. The derivation of Eq. (2) assumes that protons 
and antiprotons have the same momentum distribution 
within our acceptance as lambdas and antilambdas and 
that the decay proton or antiproton takes most of the 
momentum of the lambda or antilambda. The value of 
(p/p)/(A/A) was estimated using a coalescence model 
to relate this value to K~/K + , which is measured in 
this analysis. Values of A/p, A/p, and K~/K + were 
consistent with coalescence at 130 GeV [3,16]. A A/p 
value of 0.95 ± 0.4, based on the measured value at 
130 GeV [16], was used for the correction. The indicated 
error in this value includes the statistical and systematic 
errors in the measurement at 130 GeV (0.09 and 0.22, 
respectively), and an uncertainty associated with the 
energy dependence of this ratio. The value of K~ / ' K + 
was allowed to vary between 0.9 and 1.0 (the total 
variance allowed by our measurement). This produced 
a 1.2% correction with a systematic error contribution of 
2.2%. 

After all corrections we find the following ratios within 
our acceptance: 

(it-)/(tt + ) = 1.025 ±0.006(stat.)±0.018(syst.), 
(K-) /(K+)= 0.95 ± 0.03(stat.) ± 0.03(syst.), 
(p)/(p) = 0.73 ± 0.02(stat.) ± 0.03(syst.). 



The measured value for (p) / (p) agrees well with the 
predicted thermal model value of 0.75 at 200 GeV [17]. 
This model assumes baryon number, strangeness, and 
charge conservation and, as such, strictly applies to 4-7T 
integrated particle ratios. The prediction is based on a 
phenomenological fit to lower energy data, showing that 
the particle ratios at RHIC follow the trend observed 
at lower energies. Using the same statistical model 
calculation [18] as for our 130-GeV data [3], we estimate 
Hb/T — 0.17 ± 0.01(stat.). This estimate is consistent 
with the value of ^b/T extracted from (K + )/(K~) 
as shown in Fig. 3. For a typical chemical freeze-out 
temperature of 165 MeV, this estimate corresponds to 
Hb = 27 ± 2(stat.) MeV, which is about half of the 130- 
GeV value of fi B = 45 ± 5(stat.) MeV [3]. 
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FIG. 3. Statistical model calculation (dotted lines) of 
(K + }/(K~ ) and (p) / (p) as a function of Hb/T from Becattini 
et al. [18]. The horizontal bands show the ratios observed in 
the data (statistical errors only). The vertical shaded area 
indicates the allowed region in /xs/T. 

Our results are compared to lower energy data from 
fixed target experiments [1,2,19] and our 130-GeV mea- 
surement at the RHIC [3] in Fig. 4. Also shown are the 
results of calculations using HIJING [8] and RQMD [9] 
for 130 and 200 GeV. As was already observed at 130 
GeV [3], both HIJING and RQMD do not agree with the 
observed value of {p}/(p). 

In conclusion, the data shown in this paper provide the 
first information on the baryon density in central Au+Au 
collisions at the full RHIC energy The observed ratios 
are in agreement with expectations of statistical equilib- 
rium models. A comparison with two hadronic transport 
codes shows a significant disagreement with the ratios 
measured in this analysis. These ratios provide stringent 
constraints on current and future models dealing with 
baryon production and transport. Future analysis of p+p 
and d+Au collisions should provide information on the 
effect of multiparticlc interactions on particle production. 
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FIG. 4. (K~)/(K + ) and (p)/(p) ratios as a function of 
y?s NN for nucleus- nucleus collisions, in comparison with 
predictions from the HIJING and RQMD models. The 
percentage of the highest cross-section used for event selection 
was 12% for PHOBOS, 5% for NA49, 6% for NA44, 5% for 
AGS {K-)/{K+}, and 8% for AGS (p)/(p). Only statistical 
errors are shown. 
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